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Abstract 

We discuss theoretical uncertainties of the measurement of the W boson mass 
f£Z| at LEP2 energies, reconstructed with the help of the tandem of the Monte Carlo 

event generators KoralW and YFSWW3. Exploiting numerical results obtained 
Qh. with these programs, and the existing knowledge in the literature, we estimate 

(— i | that the theoretical uncertainty of the W mass due to electroweak corrections, as 

reconstructed at LEP2 with the help of these programs, is ~ 5MeV. Since we use 
certain idealized event selections and a simple Mp^-fitting procedure, our numerical 
exercises can be (should be) repeated for the actual ll M\y extraction methods" of the 
LEP2 measurements, using KoralW and YFSWW3 or other Monte Carlo programs. 
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In this work we would like to present our estimate of the theoretical uncertainties (TUs) 
related to electroweak corrections in the measurement of the mass of the W heavy boson 
in the LEP2 experiments. The estimate will be based on new numerical results of our 
own and on the best results available in the literature. One important reason for writing 
this paper is that the discussion of the electroweak TUs in the W mass measurement, 
including complete 0{a) electroweak (EW) corrections, is not available in the literature^]. 
On the other hand, it is becoming a burning issue, as the error on the combined LEP2 
result of the W mass measurement approaches 30 MeV, while the total TU should be 
limited to < 15 MeV. 

Since the W mass measurement has a very specific character, very different from the 
measurement of the total cross section, let us characterize it briefly. The actual way in 
which Mw is measured by the LEP2 experiments is complicated, see e.g. refs. [@-f|. In 
particular, it seems to be beyond the reach of the simple fit to a one-dimensional W 
invariant mass distributionP] (having integrated over the invariant mass of the second W). 
This is due to direct inobservability of the neutrino (from the W decay) and of most of 
the initial-state radiation (ISR) photons, loss of a fraction of the hadronic final state in 
the beam pipe and the non-trivial dependence on the invariant masses of both If' s that 
certain corrections may have. 

Let us give the reader at least a rough idea of how the W mass M w is measured by 
the LEP2 experiments. In a nutshell, this is done with the help of a two-level fitting 
procedure. At the first level, with the help of the so-called kinematic fit, an entire multi- 
momentum event, either experimental or of Monte Carlo (MC) origin, is reduced to a point 
in much fewer dimensions than the total dimension of the original set of four-momenta. 
This space consists typically of the two fitted W masses and of an auxiliary parameter 
controlling the detector energy resolution. In this way one gets 3-dimensional histograms 
with ~ 10 4 experimental events. On the other hand, one gets the analogous 3-dimensional 
histograms from a Monte Carlo simulation with ~ 10 7 events. The latter one is obtained 
typically from the combined KoralW and YFSWW3 [PHTBl programs, which we shall 
refer to as K-Y. The actual "M^-extraction" is done by fitting Mw such that the dif- 
ference between the above two 3-dimensional histograms is minimized (typically using a 
likelihood function)^. The K-Y prediction for every bin in the 3-dimensional histogram 
is of course dependent on Mw- This dependence is calculated/recalculated in the above 
fitting procedure by means of averaging the "correcting weight" P,|T3"| corresponding to 
a variation of Mw, over the entire ~ 10 7 sample of MC events stored on a computer disk. 
All complications of the experimental detector and data analysis are therefore taken into 
account without any approximation. In this way the multidifferential distribution imple- 
mented in the K-Y MC ensures a direct unbiased link between the Mw of the electroweak 
Lagrangian and the experimental LEP2 data, assuming perfect detector simulation. 

iFor instance, it is missing in the CERN Report of the 2000 LEP2 MC Workshop 0. 

2 Even more inappropriate is trying to characterize the TU of the W mass by introducing some kind 
of an "error band" in the one-dimensional W distribution - see also the discussion below. 

3 The above description of the "Mvy-extraction" tries to summarize the methods used by ALEPH, L3 
and OPAL; the DELPHI method is slightly different, see. ref. j|. 



1 



All this sounds like a strong argument to show that the theoretical uncertainties, com- 
ing from higher-order corrections and other imperfections of the theoretical calculations, 
can be studied only within the programming environment used in the actual LEP2 exper- 
iments, with the help of the K-Y MC tandem. Nevertheless, mainly because all effects 
under the following discussion are small, it makes sense to compromise and apply a "sim- 
plistic approach" based on the one-dimensional fit of a single W effective mass (integrating 
over the second one). This is what we shall do in the following. The main danger in the 
use of a fitting procedure like the one described here is that almost any physical effect 
in the W effective mass distribution may feature strong correlations as a function of the 
two effective masses, which may lead to an underestimate of the effect by a factor of 2. 
Our fitting method provides, therefore, a valuable but rough estimate of the size of the 
effects under discussion in terms of the Mw measured using LEP2 data. Consequently, if 
some effect turns out to be sizeable, that is at least 1/3 of the experimental error on Mw 
C~ 10 MeV), then it should be reanalysed within a full-scale 'Mw -extraction" procedure 
of the relevant LEP2 data analysis. In such a case, our paper can be used as a guideline 
for a more complete study to be performed by the experiments. 

Keeping all the above warnings and restrictions in mind, let us characterize more pre- 
cisely our aims and adopted methodology. During the 2000 LEP2 MC Workshop, the 
main emphasis was on the TU for the total cross section (a to t) of the M^-pair production 
process ||14|| . Since a variation of Mw is not related to the overall normalization of the 
distribution p(M 1 ,M 2 ) = da/{dM 1 dM 2 ) at .U, M, M w , but rather to the deriva- 
tives D = (d/dM i )p(M 1 ,M 2 )\ Mi=M W i the discussion of TUs on Mw is almost completely 
independent of the discussion of TUs on atot- The higher order corrections, which strongly 
influence a to t, may be completely unimportant for Mw and vice versa! In particular, it 
is inappropriate to try to translate our knowledge of TU in atot, in terms of a certain 
"error band" in the distribution da/dMi, into an error estimate of Mw - obviously it 
may easily lead to a huge overestimate of the TU of M w and to overlooking effects which 
really contribute to it. 

In the following, we consider the semileptonic process e + e~ — > udyTv^ which belongs 
to the so-called CC11 class of Feynman diagrams constituting the gauge-invariant subset 
of the 4-fermion final-state processes, see e.g. ref. for more details. We shall study 
only the leptonic W mass, i.e. the one reconstructed from the four-momenta of the \x~ 
and (in the actual experiments, the neutrino four-momentum is reconstructed from the 
constrained kinematic fit, see e.g. refs. [@-|5[]). The input parameters are the same as in the 
2000 LEP2 MC Workshop studies [ 13 1 . All the results in this paper are given for the centre- 
of-mass energy Eqms = 200 GeV and for the input W mass Mw = 80.350 GeV. The fitting 
function (FF) in all cases was taken from the semi-analytical program KorWan f|[7|,[l6|]0. 
All the results presented in the following, except the ones denoted with the label 4f , are 
for the leading-pole approximation (LPA) of the Stuart-type |17]] (the LPA a option in 
YFSWW3) applied to the above process, i.e. for the double-resonant WW production 
and decay (see ref. |L3| for more details). 

4 The relevant distribution will be available in the next release of KorWan/KoralW. 
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Figure 1: The introductory exercise, see more discussion in the text. 



In the first preparatory step, we construct a simple fitting procedure of the Mw using 
the 1-dimensional distribution of the W effective mass Mi, and we "calibrate" with the 
help of the MC data in which we switch the same effects on and off, typically the ISR and 
the non-factorizable corrections (NF) in the inclusive approximation (denoted by INF in 
the following) of the so-called screened Coulomb ansatz by Chapovsky & Khoze JTB| , just 
to see whether we get agreement in the case of the same effect in the MC data and in 
the fitting function. The other immediate profit is that we also quantify these effects as 
a shift of Mw- The results of the first exercise are shown in Fig. [|. Let us explain briefly 
the notation: Born denotes the Born-level results, ISR the ones including the 0(a 3 ) LL 
YFS exponentiation for the ISR as well as the standard Coulomb correction |f[9| , INF the 
above plus the INF correction, and Best denotes the best predictions from YFSWW3, 
i.e. all the above plus the 0{a 1 ) electroweak non-leading (NL) corrections^. 

Let us summarize observations resulting from Fig. [l|: 

• The fitted Mw exactly agrees with the input Mw in the case when the same ISR 
and INF are included both in the fitting function (FF) and the MC. 

• If one is interested only in the shift of Mw, then any of the three FFs could be used. 
In the following, in a single exercise we shall typically use one or two of them only. 

• The dependence on the fitting range is sizeable; it points out, albeit in a crude way, 
the fact discussed above, that the ultimate precise fit of Mw should always be done 
as in the LEP2 experiments, using a multidimensional fitting procedure. 



3 The 0(a 1 ) electroweak corrections for the WW production stage in YFSWW3 are based on refs. EC 
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• The size of the ISR effect is about -f MeV, that of the INF about +5 MeV, and 
the size of the NL corrections seems to be negligible, ~ 1 MeV. 

Note that there were no cuts and we used the true parton-level W invariant masses in all 
the above exercises. 

In the following exercises we shall examine the influence of various effects/corrections 
on My\r for various cuts and acceptances. Not all these effects can be included in the 
FF. Besides, only a very limited menu of cuts and acceptances can be applied in the FF. 
Therefore, our estimates of the TU will based not on absolute values of the fitted Mw but 
on relative differences of My/s corresponding to various effects. This is justified by our 
"calibration" exercise. 
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Figure 2: Effects of the ISR and the FSR on M w . 

In the second exercise, depicted in Fig. ^|, we switch on and off various orders/variants 
of the ISR and of the final-state radiation (FSR), and finally the NL correction. The FSR 
was generated by the program PHOTOS ||22|| . While the previous exercise was without 
any cuts and for the so-called BARE^ acceptance (the subscript 4ir means the full solid- 
angle coverage), we employ here the semi-realistic acceptance CAL054,,-, where all the 
photons for which the invariant mass with a final-state charged fermion was < 5 GeV 
were recombined with that fermion (also for the full solid-angle coverage). We compare 
the results for FF representing the Born and Born+ISR (no FSR) levels. 
Observations: 

• Changing the type of the ISR from 0(a 3 ) exp LL to 0(a 2 ) exp LL induces a negligible, 
< 1 MeV, effect in fitted M w . 

• The FSR effect is large for BARE 47r , ~ 60 MeV, and much smaller, ~ 7 MeV, for 
calorimetric CAL05 4vr , as expected. 
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• Switching from the single-photon (FSRi) to the double-photon (FSR2) option in 
PHOTOS results in a ~ 4 MeV change of Mw for BARE^ and no change for 
CAL05 47r . 

• The INF+NL correction is ~ 6 MeV and seems to cancel partly with the FSR, see 
CAL05 47r . 

Before we go to the next exercise, let us describe briefly the acceptances and cuts that 
were used in the MC simulations for the following calculations. 

1. We required that the polar angle of any charged final-state fermion with respect to 
the beams be 9f ch > 10°. 

2. All photons within a cone of 5° around the beams are treated as invisible, i.e. they 
were not included in the calculation of the W invariant masses. 

3. The invariant mass of a visible photon with each charged final-state fermion, Mf ch , is 
calculated, and the minimum value MJ 1 ™ is found. If MJ 1 ™ < M rec or if the photon 
energy E 1 < 1 GeV, the photon is combined with the corresponding fermion, i.e. 
the photon four-momentum is added to the fermion four-momentum and the photon 
is discarded. This is repeated for all visible photons. 

In our numerical tests we used three values of the recombination cut: 



M rec = < 5 GeV 



OGeV 
5 GeV 
25 GeV 



BARE, 

CAL05, 

CAL025. 



Let us remark that we have changed here the labelling of these recombination cuts 
from the slightly misleading bare and calo names used in Ref. |14j]. They correspond 



to our CAL05 and CAL025, respectively. This change allows us to reserve the 
BARE name for a "truly bare final fermion" setup (without any recombination). 

In the next exercise, presented in Fig. |3], we examine once again the effect of switching 
on the 4/-background corrections^ and the INF corrections, now for BARE and CAL05. 
The effect of the 4/ background is ~ 1 MeV (it is therefore negligible for the LEP2 
experiments^) and that of the INF is ~ 5 MeV. The effects of the 4/ background in the 
non-CCll channels can be larger, but they strongly depend on the applied experimental 
cuts or acceptances, so that they can be studied in detail only within the full-scale LEP2 
M\y fitting framework. 

The size of the NF effect of ~ 5 MeV requires some explanation, as the genuine NF 
effect in Mw is in fact only about ~ — 1 MeV; see refs. [f24"|-|2"7|] and the discussion of the 



6 The complete Born-level 4/ matrix element in KoralW was generated with the help of the GRACE2 
package p^j . 

7 This smallness is due to a general smallness of the 4/-background correction in the CC11 class of 
channels for LPA Q ; it may be less pronounced for a different type of LPA, such as LPAh in YFSWW3, 
for example. Our conclusion is unaffected, as it is really meant for the sum of LPA a and 4/; see also the 
discussion below. 
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Figure 3: Effects of the 4/-background, INF and NL corrections on My/. 



INF ansatz in ref. [18 (see also more discussion in the following). This effect in Fig. [3], 
understood as a difference between our ISR and INF calculations, is blown up artificially, 
because the ISR includes the ~ —6 MeV My/ shift due to the so-called "standard Coulomb 
effect" for historical reasons, although its derivation is not valid far away from the WW 
threshold. 

Another point to be explained is whether the genuine NF effect of ~ — 1 MeV obtained 
in the INF (inclusive) approximation can be increased to higher values, say 10 MeV, due 
to the LEP2 experimental cuts. In principle it can be; however, as is well known, the NF 
correction does not include (fermion) mass logarithms, and its "energy scale" , which enters 
the big logarithm owing to a cut on the photon energy, is IV and not \/s. Consequently, 
in order to get an enhancement factor of \n(Ty/ / E mauX ) ~ 10, one would need to veto 
the appearance of any photon above E max = 0.1 MeV - a very unrealistic experimental 
selection indeed. On the contrary, in the actual LEP2 experiments photons with energy 
< 2 GeV are not disturbed, directly or indirectly, by the experimental event selection. 
This is why any strong enhancement of the NF effect with respect to its "inclusive" 
treatment (INF) must be just absent. 

In the Fig. |4], we examine the difference in My/ fitted to the W mass distribution 
obtained from YFSWW3 and from Racoon WW |28, 29fl. The distributions used for the 
My/ fits are exactly the same as those that were used for the plots in the CERN Report 
of the 2000 LEP2 MC Workshop, see [rj4 ]. The statistical error is taken into account 
in the fits and propagates into the fitted W mass. It is merely < 1 MeV. We use two 
fitting functions, one in which the ISR is included (with the incomplete NL but with the 
Coulomb effect) and another one in which the INF (Chapovsky & Khoze) is also included. 
It is done for two kinds of calorimetric acceptances: the CAL05 and CAL025 described 
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Figure 4: Comparison of YFSWW3 and RacoonWW. 

above. 

Observations concerning the results shown in Fig. |j: 

• The comparison of YFSWW3 with RacoonWW is very interesting because the two 
calculations differ in almost every aspect of the implementation of the ISR, FSR, 
NL and NF corrections. 

• It is quite striking that the results of YFSWW3 and RacoonWW differ, in terms of 
the fitted mass, by only < 3 MeV, slightly more for CAL05 than for CAL025. 

• The difference between YFSWW3 and RacoonWW is definitely smaller than the 
size of the INF correction, roughly by a factor of 2 (the INF is of order 3-5 MeV 
for these two acceptances). 

The most important result of the comparison between YFSWW3 and RacoonWW is that 
it reconfirms the smallness of the NF corrections in the W mass. Its size is well below 
the 10 MeV precision target of the TU for the measurement of at LEP2. It would be 
interesting to repeat the above exercise for the true LEP2 acceptance, using the full-scale 
fitting procedure. In our opinion the above difference between YFSWW3 and RacoonWW 
in terms of the fitted Mw cannot be attributed to some dominant source. Most likely 
it consists of several contributions and the leading candidates are ordinary factorizable 
QED corrections and/or some purely technical/numerical problems. 

In order to gain better understanding of the above numerical results it is worth while 
to estimate them semi-quantitatively, in terms of some "scale parameters" representing 
various QED or EW corrections to Mw- We shall do it in the following. In addition, we 
shall discuss certain effects not included in K-Y or RacoonWW. The corrections to Mw 
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AM W 


Error Type 


Scale Param. AMw = T x e 


Numerical cross-check 


AM W 


WW production 


ISR 0(a 4 Lj) 
ISR 0{a 2 L e ) 
ISR 0(a 2 ) pairs 


e~ v -^P^( a fLt ~ 5- 1(T 6 
e~ r ^(«)2 L 2 _ 4 . 1Q -4 


[0(a 3 L 3 ) - 0(a 2 L 2 e )] Koralw 
KorWan 
KorWan 


< 1 MeV 

< 1 MeV 

< 1 MeV 


W decay 


FSR 0(a) miss . 


e ~ 0.2 f ^2 In ~ 10" 3 


Basic tests of PHOTOS 


~ 2 MeV 


FSR 0(a 2 ) mtss . 


e~|ff^21n^) ~ 10~ 5 
2 \ 8 7r pt y 


On/off 2 7 in PHOTOS 


< 1 MeV 


Non-factorizable QED interferences (between production and 2 decays) 


0((X )miss. 


e ~o.i(5 (1 -" >a )~io- 4 


Chapovsky & Khoze 


< 2 MeV 


0(a 2 ) 




None 


< 1 MeV 



Table 1: Estimation of the missing effects in the K-Y MC tandem. 



are generally of the type ~ or ~ (j^-) e, where £ is a small parameter of 
the perturbative expansion. We divide them into three types: 

Case (a): A mildly mass-dependent correction to the W mass distribution p(M) = 4?7 ~ 

\BW{M)\ 2 x f(M 2 ), which leads to 5M W ~ |r 2 w dln /|f 2) \ M=Mw , where SW(M) 
denotes the Breit-Wigner resonance function and f{M 2 ) is a mild function in the 
vicinity of the resonance (in the semi-quantitative discussion, we usually take M = 
(M\ + M2)/2). The most trivial example is the kinematic factor f(M 2 ) = (3m = 
(1 - AM 2 /s) l ' 2 } yielding 5 kin M w ~ _r w ^^, where /3 W = /3 m |m=m w - It is 
not visible in our fits (always taken into account in the FF); however, the f3w~ 
factor gets modified by the ISR, giving rise to 5M W ~ 5ki n M w x 2^L e ~ —6 MeV 
(L e = 21n^-)Q This effect is responsible for most of the M\y shift when switching 
from the Born to the ISR in Fig. |l| It vanishes at high energies. The response of 
SMw to a more general variation: f(x) — > f(x) + efi(x), where efi is due to the 
higher-order ISR correction, is in general negligible, < 1 MeV. It can be estimated 

using 5M W ~ ^\^w dXn ~ \m=m w — £ (here, we exploit the fact that fx has 
a derivative of 0(1) as a function of M 2 /s). For instance, the missing 0(a 2 ) NLL 
ISR is proportional to e ~ a 2 L e ~ 10~ 3 , giving rise to 5Mw ~ 10 -3 MeV. 

Case (b): This is the case of the QED effects in the decays, the so-called FSR. In this 
case the mass distribution gets distorted according to -^^(M 2 ) ~ J dz Jfsr(^ — 
z)^ FSR - 1 ^(zM 2 ), where -? FSR ~ (f)ln(M^/mJ) ~ 0.03 for the BARE and 

8 It is obtained from the approximate evaluation of the derivative of the ISR convolution: 
B& [f dz^li SR (l - z)™*- 1 ] |m=m wi where 1ISR = %L e . 
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Ifsr — (^)^ n ( M w/ m CALo) — °- 01 for CALO-type acceptance (with m CALO = 
5 GeV). The mass shift 5M W ~ Tyye, e ~ — | Jfsr — —0.012, is accounted for in 
the complete 0(a) calculation; see also refs. fpO|, plfl • In the case of PHOTOS the 
missing 0(a) is related mainly to high-pp photons, and from tests of this program 
listed in ref. JEJ one can conclude that it corresponds to ~ 0.2 x e, The 
missing 0(a 2 ) FSR effect we estimate as follows: A5M W ~ T w \e 2 <C 1 MeV. 



Case (c): The influence of the NF QED interferences on the W mass is character- 
ized by the correction function 5nf, which is strongly dependent on My/ in the 

vicinity of the resonance: ~ \BW(M)\ 2 f(M) l + a5 NF (- 



m w t 



The re- 



sulting Mjy-shift is 5M W ~ \^w a ^^M^'\M=M w - The simple INF formula for 



W« dM \M=M W 

S NF of ref. fTJ, leads to 5M W ~ -IVf 1 MeV, in perfect agree- 



ment with the complete calculations of the NF corrections, see refs. [p^-p7|. If 
0{a 1 ) NF is accounted for, then we estimate the missing 0(a 2 ) NF contribution at 

A5M W ~ Y w \ ^^)iy « 1 MeV. 

The above discussion confirms that all our numerical results are consistent with expec- 
tations based on the "scale parameters" analysis and semi-quantitative calculations, and 
provides some estimates of the effects not accounted for in our MC programs. It is also 
summarized in Table [I] together with the relevant numerical estimates. 

Let us finally discuss a question of the TU due to the so-called "ambiguity of a defi- 
nition of the LPA" . In YFSWW3, we implemented two different definitions of the LPA, 



called LPA a and LPA?, fL3| . Differences in My/ reconstructed from the results obtained 
in these two options can give us a hint of the TU due to the LPA. Actually, we need to 
check only the variation of My/ caused by the NL corrections. This is because KoralW 
implements the full 4/-process at the so-called ISR level; hence, the ambiguity due the 
LPA is reduced from 0(j^-) to j^-) and is located only in the NL part. We per- 
formed numerical tests of the dependence of NL on the choice of the LPA with the help 
of YFSWW3, finding the variation of the My/ induced by the change from LPA a to LPAf, 



to be < 1 MeV (the ISR and NL parts are always defined as in refs. || 13]). Another 
uncertainty in the LPA is due to the missing higher orders in the NL part. This can 
be estimated by comparing the predictions of the so-called schemes (A) and (B) ]13[ in 
YFSWW3. These two schemes account for some higher-order effects by the use of the 
effective couplings in two different ways - in fact, the scheme (B) follows the prescription 
employed in RacoonWW. We have checked that the change from the scheme (A) to the 
scheme (B) results in the fitted W mass shift of < 1 MeV (as expected, the results of the 
latter scheme are slightly closer to the ones of RacoonWW). Consequently, we attribute 
AM W = 1 MeV to the TU of M w due to the LPA. 

In the above, we have considered only the leptonic W mass coming from the W~ — > 
yU~z^ decay. For the other leptonic decays, the results should be similar when one applies 



the calorimetric-type acceptance - as was shown in ref. [11]. For the hadronic W masses at 



the parton level, the results would be analogous to the ones presented here with the FSR 
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switched ofiQ. More realistic estimates of the TU for the hadronic W mass would require 
taking into account the QCD effects, hadronization, jet definitions, etc. This should be 
done in the full-scale experimental data analysis, which is beyond the scope of this paper. 

From the above numerical exercises and the accompanying discussion, we come to the 
following conclusions: 

• The electroweak theoretical uncertainty in Mw of the K-Y MC tandem at LEP2 
energies is ~ 5 MeV. 

• The above conclusion is strengthened by the smallness of the differences between 
YFSWW3 and RacoonWW, which we attribute to the standard factorizable correc- 
tions (ISR, FSR, etc.) and purely technical/numerical effects. 

• In the above estimate we included a "safety factor" of 2, corresponding to the fact 
that our fits of Myy were done for 1-dimensional effective W mass distributions. In 
order to eliminate it, our analysis should be repeated for the realistic measurements 
of the LEP2 experiments. 
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